WO 2005/036758 ™ *» ™* « "Vff/Ht^ * ^" " 



KR20047002630 

1 - 



METHOD FOR ENCODING LOW-DENSITY PARITY CHECK CODE 
BACKGROUND OF THE INVENTION 

. L Field of the Invention 

The present invention relates generally to a method for encoding data, 
and in particular, to a method for encoding a low-density parity check (LDPC) 
code. 

2. Description of the Related Art 

In general, a communication system encodes transmission data before 
transmission to increase stability of transmission, and prevents excessive 
retransmissions to increase transmission efficiency. For coding the transmission 
data, the mobile communication system uses Convolutional Coding, Turbo 
Coding, and Quasi Complementary Turbo Coding (QCTC). The use of the 
coding schemes stated above contributes to an increase in stability of data 
transmission and transmission efficiency 

Recently, wireless communication systems are evolving into advanced 
wireless communication systems capable of transmitting data at very high speed. 
The advanced wireless communication system desires to transmit data at higher 
speed. Accordingly, there is a demand for an advanced coding scheme capable of 
obtaining higher efficiency than that of the current coding schemes stated above. 

Low-density parity check (LDPC) coding is provided as a new coding 
scheme to meet the demand. A detailed description of the low-density parity 
check code will now be described herein below. The low-density parity check 
code was first proposed by Gallager in the early 1960s, and reviewed by 
MacKay in the late 1990s. The low-density parity check code reviewed by 
MacKay is based on sum-product algorithm. Since the use of belief propagation 
decoding, die low-density parity check code has started to attract public attention 
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as a code capable of showing excellent performance approximating the Shannon 
capacity limit. 

Thereafter, Richardson and Chung proposed a density evolution 
5 technique for tracing a variation according to iteration in probability distribution 
of messages generated and updated during decoding on a factor graph 
constituting a low-density parity check code. For the density evolution technique 
and infinite iteration on a cycle-free factor graph, Richardson and Chung 
invented a channel parameter (or threshold) capable of enabling error probability 
10 to converge into '0\ That is, Richardson and Chung proposed degree distribution 
capable of maximizing channel parameters of variable nodes and check nodes on 
the factor graph. In addition, Richardson and Chung theoretically showed that 
such a case can be applied even to an LDPC code with a finite length in which 
there are cycles. 

15 

In addition, Richardson and Chung showed that theoretical channel 
capacity of an irregular LDPC code can approximate the Shannon capacity limit 
up to only 0.0045 dB using the density evolution technique. In particular, Flarion 
Co., leading realization of design and hardware (H/W) of the LDPC code, has 
20 proposed a multi-edge type vector LDPC code capable of realizing a parallel 
decoder having a frame error rate lower than that of a turbo code even for an 
LDPC code with a short length. 

The LDPC code is treated as a powerful alternative to the turbo code in 
25 the next generation mobile communication system. This is because of parallel 
structure and low complexity of the LDPC code for decoder realization, and low 
error floor and good frame error rate in terms of performance. Therefore, it is 
expected that the future researches will provide LDPC codes having better 
characteristics. 

30 

However, in realization, the current LDPC code is more complex than 
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the turbo code in encoding process, and requires a structure of an optimized code 
capable of providing better performance than the turbo code at a short frame size. 
Although active researches have been made to solve this problem, there has been 
proposed no scheme capable of encoding an optimized LDPC code. 

SUMMARY OF THE INVENTION 

It is, therefore, an object of the present invention to provide an LDPC 
coding method with a simple encoding process. 

It is another object of the present invention to provide an LDPC coding 
method having improved performance at a short frame size. 

To achieve the above and other objects, there is provided a method for 
generating a low-density parity check code consisting of an information-part 
matrix and a parity-part matrix. The method comprises the steps of changing the 
information-part matrix to an array code structure, and allocating a degree 
sequence to each of submatrix columns; extending the parity-part matrix such 
that an offset value between diagonal lines has a predetermined value in a 
generalized dual-diagonal matrix which is the parity-part matrix; lifting the 
generalized dual-diagonal matrix; determining an offset value for cyclic column 
shift for each submatrix of the lifted generalized dual-diagonal matrix; and 
performing an encoding process for determining a parity symbol corresponding 
to a column of the parity-part matrix. 

Preferably, the degree sequence consists of 

D = {15, 15, 15, 5, 5, 5, 4, 3, 3, 3, 3, 3, 3, 3, 3} 

, and the offset value between diagonal 
lines is relatively prime to the number of columns in the generalized dual- 
diagonal matrix. 
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Preferably, a difference between a sum of offset values for cyclic row 
shift of a submatrix on a diagonal line in a generalized dual-diagonal matrix 
which is the parity-part matrix and a sum of offset values for cyclic row shift of 
a submatrix on an offset diagonal line is not 0. 

5 

Preferably, the encoding process comprises the process of (a) 
determining a parity symbol of a first row in a submatrix with a submatrix 
column index 0 on a diagonal line of the parity-part matrix; (b) setting a row 
index in a submatrix of a parity symbol being identical to the determined parity 

1 0 symbol in column index in a submatrix in a submatrix on an offset diagonal line 
having the same submatrix column index as the submatrix column index of the 
set parity symbol; (c) determining a parity symbol having the same row index in 
the set submatrix in a submatrix on a diagonal line having the same submatrix 
row index as the submatrix row index of the submatrix on the offset diagonal 

15 line; and (d) repeatedly performing the steps .(b) and (c) until generation of the 
parity matrix is completed. 

Preferably, in step (a), the parity symbol is determined by a sum of 
information symbols of the information-part matrix existing in the same row as a 
20 row index in the submatrix whose parity symbols are determined. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The above and other objects, features and advantages of the present 
25 invention will become more apparent from the following detailed description 
when taken in conjunction with the accompanying drawings in which: 

FIG 1 is a diagram illustrating a parity check matrix for a general (p,r) 
array code; 

FIG. 2 is a diagram illustrating an example of a matrix H d in the case 
30 where the maximum number of Is existing in a particular column is defined as 
d v , the number of such submatrix columns is defined as n v , and the number of Is 
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existing in the remaining submatrix columns is always 3; 

FIG. 3 is a diagram illustrating a factor graph structure of an irregular 
repeat accumulate code; 

FIG. 4 is a diagram illustrating a matrix of a low-density parity check 
code having an irregular repeat accumulate code; 

FIG. 5 is a diagram illustrating a parity matrix in which an offset value 
T of a dual-diagonal matrix is lifted to a particular value; 

FIG 6 is a diagram illustrating a process of sequentially calculating 
values of p 0 , p,, f , p r . 2f , - in the encoding process of FIG. 5; 

FIG 7 is a diagram illustrating an example of an information part H ( , of 
a parity check matrix generated in a method described in section A, wherein a 
maximum variable node degree is 1 5; 

FIG 8 is a diagram illustrating a 12x12 matrix obtained by matrix-lifting 
a basic 4x4 matrix in a method of replacing each element in a 4x4 matrix with a 
3x3 identity matrix or a 3x3 0 matrix; 

FIG. 9 is a diagram illustrating a parity matrix H p constructed through 
matrix lifting by a pxp cyclic permutation submatrix; 

FIG. 10 is a flowchart illustrating a method for generating a matrix of a 
parity part according to a preferred embodiment of the present invention; 

FIG. 11 is a diagram illustrating a parity matrix of a generalized dual- 
diagonal matrix lifted for r=l 5, f=7 and p=89. 

FIG. 12 is a flowchart illustrating an iterative belief propagation 
decoding process for proving efficiency of the present invention; 

FIG. 13 A is a diagram illustrating an example of a matrix of an 
information part H d for n=870 and p=29; 

FIG 13B is a diagram illustrating an example of a matrix of a parity pari 
H p forn=870andp=29; 

FIG. 13C is a simulation result graph illustrating a comparison between a 
low-density parity check code and a turbo code; 

FIG. 14A is a diagram illustrating an example of a matrix of an 
information part H u for n=1590 and p=53; 
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FIG 14B is a diagram illustrating an example of a matrix of a parity part 
Hp for n=1590 and p=53; 

FIG. 14C is a simulation result graph illustrating a comparison between a 
low-density parity check code and a turbo code; 
5 FIG. 15A is a diagram illustrating an example of a matrix of an 

information part H d for n=3090 and p=103; 

FIG. 1 5B is a diagram illustrating an example of a matrix of a parity part 
H p for n=3090 and p=103; 

FIG. 15C is a simulation result graph illustrating a comparison between a 
1 0 low-density parity check code and a turbo code; 

FIG 16A is a diagram illustrating an example of a matrix of an 
information part H d for n=7710 and p=257; 

FIG. 16B is a diagram illustrating an example of a matrix of a parity part 
H p for n=7710 and p=257; 
1 5 FIG 16C is a simulation result graph illustrating a comparison between a 

low-density parity check code and a turbo code; and 

FIG 1 6D is a simulation result graph for a low-density parity check code 
according to a variation in the number of iterations. 

20 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

A preferred embodiment of the present invention will now be described 
in detail with reference to the annexed drawings. In the drawings, the same or 
similar elements are denoted by the same reference numerals even though they 
25 are depicted in different drawings. In the following description, a detailed 
description of known functions and configurations incorporated herein has been 
omitted for conciseness. 

In the following description, the present invention proposes a new 
30 LDPC code capable of showing simple coding and good performance. To this 
end, the present invention defines a new parity check matrix. That is, the present 
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invention proposes a new LDPC code that uses a new parity cheek matrix. In 
addition, the present invention shows that the new LDPC code can be simply 
encoded by linear calculation, and proposes a coding method for the same. 
Finally, the present invention shows that due to belief propagation decoding, the 
new LDPC code has better decoding performance than that of the turbo code 
used in the CDMA2000 lx EV-DV standard. 

For convenience, a description the present invention will be made as to 
performance and realization of an LDPC code with a code rate of 1/2. However, 
the code rate can be lifted within the scope of the present invention. 

1. Parity Check Matrix Desip n 

In this section, the present invention defines a new parity check matrix 
generated by grafting and applying parity check matrix structures that define an 
array code and an irregular repeat accumulate (IRA) code. Also, a description 
will be made of a method for generating a matrix having a larger size while 
maintaining a characteristic of a basic matrix generated in this manner. 

A. Array Code Structure 

Generally, a parity check matrix for a (p,r) array code is defined as FIG. 
1. FIG. 1 is a diagram illustrating a parity check matrix for a general (p,r) array 
code. With reference FIG. 1, a description will now be made of a parity check 
mahix for a general (p,r) array code. 

In FIG. 1, p denotes a prime number, and a* denotes a pxp cyclic 
permutation matrix obtained by cyclic-shifting respective rows of a pxp identity 
matrix I which is a square matrix having a size of p, by j. A row and a column 
each composed of a set of the cr* are called a submatrix row and a submahix 
column, respectively. A column and a row of a parity check matrix defining an 
array code uniformly have p Is and r Is, respectively. When the p is reasonably 
large, a ratio of Is in the matrix is reduced, making a structure of a low-density 
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parity check code. Such an array code does not have a cycle structure with a 
length of 4. That is, if elements belonging to 4 partial matrixes a 13 , a ,b , a* 3 and 
<yi b (i^j) forming a square, or cycle 4, in a parity check matrix have a cycle-4 
structure, the following relation should be satisfied. 

5 

ia-ib=ja-jb Q (i-j) = (a-b) (mod p) 



In Equation (1), it should be satisfied that a = b. However, because 'a' 
and 'b' exiting in different rows always have different values, Equation (1) 
10 cannot be satisfied. Therefore, an array code having the parity check matrix 
described above has not a cycle-4 structure. 

The present invention generates a new matrix structure by modifying the 
parity check matrix structure of the array code described above. A description 

1 5 will now be made of a process of generating a new matrix structure to be used in 
the present invention. The object of the modification is to obtain an irregular 
structure in which distribution of 1 existing in each column of a matrix to be 
generated is irregular, based on an array code structure not having a cycle-4 
structure. This process should be performed such that distribution of 1 existing 

20 in each row is relatively uniform, i.e., only a total of 2 types exist. A description 
will now be made of a method for generating a new matrix structure used in the 
present invention. 

(1) A j th submatrix column of a parity check matrix forming an array 
25 code is composed of submatrixes represented by 



//.=[/ a' * 2J - <j {r - l)J ) r 



(2) 
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(2) Depending on a predetermined degree distribution, a degree 
sequence of each submatrix column is defined as 



W/ 



iii 



In Equation (3), dj denotes a column degree corresponding to a j 
submatrix column, and V denotes the total number of submatrix columns. In 
addition, the V is set such that it is identical to the total number, d v , of Is 
existing in a particular column. 



(3) Depending on the degree sequence D defined in (2), a j ,h submatrix 
column is modified as 



/4=[0 0 - a** a''""" - ***** 0 - Of, t, = /,_, + d,. x (mod S ) 



In Equation (4), tj denotes a submatrix row number from which a non-0 
submatrix starts in a j th submatrix column Hj. 

(4) A matrix H d to be generated is defined by Equation (5) in which 
submatrix columns are replaced with column vectors. 

H d =[// 0 H y - Hj ••• //,.,] 

When a particular submatrix column of a matrix is defined according to 
the process described above, a particular column existing in each submatrix 
column always has as many Is as only a degree allocated to the corresponding 
submatrix column. A matrix composed of the submatrix columns will minimize 
submatrix overlap between matrixes. Therefore, the number of Is existing in a 
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particular row of the entire matrix will also be minimized. That is, it can be 
easily understood that the number of Is in one row in the entire matrix will 
always be 2 or 1. Because the generated matrix H d has a structure obtained by 
removing some submatrixes from the array code structure, it always does not 
5 have a cycle-4 structure like the array code. 

FIG. 2 is a diagram illustrating an example of a matrix H<i in the case 
where the maximum number of Is existing in a particular column is defined as 
d v , the number of such submatrix columns is defined as n v , and the number of Is 
10 existing in the remaining submatrix columns is always 3. 

In the present invention, an information-part matrix H d 220 of a parity 
check matrix H defined based on the array code structure described in 
connection with FIG. 2 will be used as a submatrix of a parity check matrix H of 
15 a low-density parity check code to be defined in the present invention. A 
description will now be made of an irregular repeat accumulate (IRA) code for 
generation of a submatrix forming a parity part of the parity check matrix H. 

B. Generalized Dual-Diagonal Matrix 
20 It is known that an irregular repeat accumulate (IRA) code has a simple 

encoder structure and has relatively high performance by a message passing 
decoder. Such an irregular repeat accumulate code can be regarded as one type 
of the low-density parity check mode. 

25 FIG. 3 is a diagram illustrating a factor graph structure of an irregular 

repeat accumulate code. With reference to FIG. 3, a description will be made of a 
factor graph structure of an irregular repeat accumulate code. The factor graph of 
the irregular repeat accumulate code of FIG 3 is a systematic version, in which 
all parity nodes corresponding to parity symbols of a codeword except only one 

30 parity node have a degree of 2. That is, each of parity nodes 301, 302, 303, 

304, 305, 306 and 307 has 2 edges connected thereto. Such a parity node 
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structure is advantageous in that generation of parity symbols based on given 
information symbols is simply achieved by linear calculation. In addition, as can 
be understood from FIG. 3, information nodes of the irregular repeat accumulate 
code can have various degree distributions. Values of the information nodes are 
connected to check nodes through a random permutation 320. 

A matrix construction of the irregular repeat accumulate code having the 
factor graph of FIG. 3 can be illustrated in FIG. 4. FIG. 4 is a diagram illustrating 
a matrix of a low-density parity check code having an irregular repeat 
accumulate code. 

In the matrix of a low-density parity check code illustrated in FIG. 4, 
diagonal lines in a parity part represent T, and an offset «f (401) between the 
diagonal lines have a value of 1 in a general irregular repeat accumulate code. 
Therefore, a submatrix corresponding to the parity part of a parity check matrix 
defining an irregular repeat accumulate code has a construction of a dual- 
diagonal matrix. The number of Is existing in each column of an information 
part in the parity check matrix has irregular distribution. The distribution is 
defined by a random permutation unit of FIG. 4. 

The present invention considers only the parity part of the irregular 
repeat accumulate code. Therefore, the present invention generalizes only the 
parity part and proposes a new type of matrix. 

In the case of the general irregular repeat accumulate code of FIG. 4, an 
offset T between diagonal lines has a value of 1 (f=l). In the present invention, 
if the offset value T is lifted so that it has a particular value, a dual-diagonal 
matrix of the parity part can be generalized as illustrated in FIG. 5. In FIG. 5, a 
first diagonal line 501 forms a diagonal line from a position in the first column 
and first row to a position in the last column and last row in the matrix. A first 
partial diagonal line 502a of a second diagonal line has a value of '0' in a first 
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position column-shifted by the offset value, and has a value of T in the 
remaining positions up to the last point thereof. A second partial diagonal line 
502b of the second diagonal line forms a diagonal line with a value of c V from a 
first column in the next row. 

5 

Herein, FIG. 5 is a diagram illustrating a parity matrix in which an offset 
value of a dual-diagonal matrix is lifted to a particular value. Although the offset 
value is changed as illustrated in FIG. 5, there always exists a column having 
only one c 1 \ It should be noted that if the number of Is existing in all columns is 

10 2, a rank of the dual-diagonal matrix becomes smaller than the number of rows 
in the matrix. Therefore, the first row of the second diagonal line should have c 0' 
rather than T as illustrated in FIG 5. In a case of an rxr matrix, i.e., a square 
matrix with a matrix size c r 5 illustrated in FIG. 5, if an offset value T of a dual- 
diagonal matrix and a size c r' of the matrix have no common factor, it is possible 

15 to generate all of r parity symbols by simply performing addition calculation on 
a set of given information symbols r times. In order to prove this, consideration 
will be taken into the following theorems. 

Theorem 1: In an Abelian group AG r ={0,l,2,-,r~l}, all elements are 
20 generated by performing addition on a particular element T which is relatively 
prime to V and is not '0', less than V times. In order to prove this, if all 
elements in the AG r cannot be generated by performing addition on the particular 
element T which is relatively prime to V and is not '0', less than 'r 5 times, 
there is it' less than 'r', which satisfies a condition of Equation (6). 

25 

kf= 0 (mod i) f ~ 



30 



However, Equation (6) is against the assumption that T is relatively 
prime to c r\ That is, a minimum value of the 'k 5 satisfying Equation (6) should 
always be 6 r\ Therefore, all elements in the AG r can be generated by performing 
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addition calculation on the element T r times. 

The dual-diagonal matrix illustrated in FIG. 5 is a matrix corresponding 
to a parity part of a parity check matrix. Therefore, as to an i* column in the 
matiix of FIG. 5, it is possible to always obtain information shown in Equation 
(7) from the i* column in a matrix corresponding to the remaining information 
part of the parity check matrix. 



In Equation (7), d; denotes an i' h information symbol., and a { denotes the 
number of Is existing in an i ,h column of a matrix corresponding to an 
information part of a parity check matrix. In addition, I denotes addition 
calculation in GF(2). The GF(2) refers to Galois Field over 2, and means a finite 
field defined in modulo 2. The term "finite field" refers to a set, in which the 
number of elements is finite and the elements are closed for addition and 
multiplication over modulo 2, there are an identity element and an inverse 
element for addition, there are an identity element for multiplication and an 
inverse element for a non-0 element, and commutative low, associative law and 
distributive law are satisfied for the addition and multiplication. Therefore, the 
GF(2) means a set of {0,1} satisfying such a property. As a result, a first parity 
symbol that can be obtained by a parity check equation corresponding to a first 
column in the matiix of FIG. 5 is defined as 



In the matrix construction of FIG. 5, a relation between a column index 
and a row index for particular Is existing in a diagonal line and an offset 
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diagonal line is defined as 



If two formulas shown in Equation (9) are used, p r . f can be calculated 
using p 0 in accordance with Equation (10). 



Pr- f=n + K-r 



(10) 



In addition, because a parity check matrix has a construction of a dual- 
diagonal matrix as illustrated in FIG. 5, Equation (10) can be modified into 



Pr-U ~ Pr-(l-n' + V r-V ' = L-.r-l ^ 



111 this way, a particular parity symbol p r . if can be obtained by Equation 
(11), and based on Theorem 1, if V is relatively prime to T, as £ i' increases, a 
value of (r-if) (mod r) has all value of 1 to r-1 only once. Accordingly, if V is 
relatively prime to T, values of all parity symbols can be calculated through the 
foregoing process. Therefore, in the foregoing method, the parity check matrix 
having the parity part illustrated in FIG. 5 can always undergo simple encoding 
by linear calculation. FIG. 6 is a diagram illustrating a process of sequentially 
calculating values of p 0 , p r . f , p r . 2f , ••• in an encoding process. 

Referring to FIG. 6, a connection is made from a position value of a 0 th 
column and a 0 th row to a position value of the 0th column and an (r-f) ,h row as 
represented by reference numeral 601. It can be noted that the position value of 
the 0 th column and the (r-f) th row is connected to a position value of an (r-f)" 1 
column and the (r-f) th row as represented by reference numeral 602. In this way, 
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ls are located only in a first diagonal line 610, a first partial diagonal line 620a 
of a second diagonal line, and a second partial diagonal line 620b of the second 
diagonal line. Because only one value of 1 exists in a first column of the 
diagonal lines, a second-order function is calculated by continuously finding a 
value of 1 in the forgoing manner. As a result, all values can be found. 

The generalized dual-diagonal matrix construction described above will 
be used as a submatrix H p of a parity part in a parity check matrix H of a low- 
density parity check code to be defined in the present invention. A description 
will now be made of a method for generating a parity check matrix of a low- 
density parity check code to be defined in the present invention. 

C. Parity Che ck Matrix Construction 

In this section, a description will be made of a method for generating a 
new parity check matrix H according to the present invention using a 
generalized dual-diagonal matrix construction described in section B, and also 
made of a construction thereof. A parity check matrix described below defines 
matrix H„ defined based on the array code structure described in section A, 
information part of a parity check matrix H. Further, the parity check matrix 
described below defines the normalized dual-diagonal matrix H p described in 
section B as a parity part of the information-part matrix H. Therefore, a parity 
check matrix H defining a low-density parity check code to be designed in the 
present invention has a construction represented by 

H = [H d | H„] 

where H d : parity check sub-matrix for information part 
parity check sub-matrix for parity part 

(12) 

In Equation (12), respective submatrix columns in an information-part 
matrix H d of the parity check matrix are defined by a degree sequence D based 



a 

as an 
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on a predefined degree distribution, and a unit submatrix a ,J constituting each 
submatrix column is a cyclic permutation matrix obtained by column-shifting a 
pxp identity matrix with a size of £ p' by ij which is a product of a submatrix 
column index 'i' and a submatrix row index c j\ Here, 'p' is always a prime 
5 number. In addition, H p lifts Is and Os existing in a matrix to a pxp matrix and a 
0 matrix, respectively, in order to match its construction to a construction of H d . 

A description will now be made of processes of generating a low-density 
parity check code according to an embodiment of the present invention. 

10 

(I) Degree Distribution and H,i Construction 

Richardson et al. showed that when message passing decoding is 
achieved in a factor graph defining a low-density parity check code, it is possible 
to trace a variation in messages by a variable & check node message update 

1 5 process and an iterative processing process by a sum-product algorithm through 
a variation in probabilistic distribution. Further, Richardson et al. showed that 
there is a chamiel parameter threshold capable of determining whether an 
average error probability of a low-density parity check code defined with a 
factor graph having a degree distribution through such a density evolution 

20 technique converges into c 0\ In addition, Richardson et al. showed that for a 
channel parameter lower than the channel parameter threshold, when infinite 
iteration is assumed in a decoding process, probability of bit error always can 
converge into c 0\ Therefore, such a density evolution technique can be used as a 
design tool capable of optimizing degree distribution of a low-density parity 

25 check code capable of improving a threshold for a particular channel 
environment. Thus, many researches have been made to calculate a degree 
distribution optimized for a low-density parity check code and a threshold at that 
time. 

30 The present invention approximates a predefined optimized variable 

node degree distribution to a variable node degree distribution of a parity check 
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matxix H to be defined in the present invention. In this case, variable nodes with 
a degree of 2 are all allocated to respective columns in a parity part H,, of the TJ, 
and variable nodes with a degree higher than 2 are allocated to respective 
columns in an information part H d of the H. In addition, from the array code 
structure described in section A and the submatrix construction of the 
generalized dual-diagonal matrix described in section B, a check node degree of 
the parity check matrix H always has only one or two minimized types in terms 
of the number of all given Is. 



Now, a degree sequence for defining an information part H d of a parity 
check matrix to be defined in the present invention will be determined. If a 
maximum variable node degree in a factor graph defining a low-density parity 
check code is defined as d v , distribution of respective edges existing in a variable 
node can be expressed with a polynomial represented by 

(13) 

In Equation (13), X ( denotes a ratio of edges existing in a variable node 
with a degree of e i\ A ratio Cj of edges existing in a variable node with a degree 
of 'j' in a factor graph for a given X(x) can be expressed as 



c j = 



te2 



(14) 



Considering implementation in Equation (14), a construction of a parity 
check matrix will be described for a maximum variable node degree = 15. It is 
known that when a value of a maximum variable node degree in a factor graph 
defining a low-density parity check code is 15, an optimal degree distribution for 
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a code rate of 1/2 in a binary input additive white Gaussian noise (Bi-AWGN) 
channel environment is given as shown in Equation (15) by the density evolution 
technique. 

X 2 =0.23802, X 3 =0.20997, X 4 =0.03492, X 5 =0.12015, 
X 7 =0.01587, X X4 =0.00480, X l5 =0.37627 



In an additive white Gaussian noise (AWGN) channel for allowing error 
probability of variable nodes to converge into c 0\ a maximum channel 
parameter noise variance a* is 0.9622, and this becomes 0.3348 dB in ternis of 

1 0 Eb/No. That is, a low-density parity check code defined by a factor graph having 
the foregoing degree distribution shows performance approximating the 
Shannon capacity limit by 0.3348 dB by belief propagation decoding for which 
infinite block size and infinite iteration are assumed. In a factor graph having the 
foregoing degree distribution, ratios of variable nodes having respective degrees 

1 5 are given as 

Cj = 0.47709, c 3 = 0.28058, c K = 0.034997, c 5 = 0.096332 
c, = 0.0090386, c H = 0.001374, q 5 = 0.100560 

(16) 



In the parity check matrix to be defined in the present invention, variable 
20 nodes with a degree of 2 should always be allocated to a parity part H p , and only 
other variable nodes except the variable nodes with a degree of 2 should always 
be allocated to an information part Hj. In addition, all columns in each 
submatrix column of the information part H d have the same column weight. 
Therefore, ratios of variable nodes having respective degrees in a factor graph of 
25 the parity check matrix H are approximated as shown by Equation (17). 
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c 2 = 1/2 « 0.5, c 3 = 8/30 = 0.26667, c 4 = 1/30 - 0.03333, c 5 = 3/30 - 0.1 , 
q 4 = 0, c l5 = 3/30 = 0.1 



With use of the variable node ratios shown in Equation (17), a degree 
distribution polynomial A,(x) is calculated gain, and a threshold a* in AWGN of a 
low-density parity check code defined by the corresponding degree distribution 
is calculated in accordance with Equation (18). 



a =0.9352 (Eb/No)' = 0.5819 dB 



A value determined by Equation (18) is a value corresponding to a low- 
density parity check code with a code rate of 1/2. This value was derived using 
the density evolution technique based on a Gaussian approximation value as an 
analysis technique, and if estimated probability of bit error is less than 1 0" 6 , the 
corresponding value is calculated considering that there is no error. A degree 
sequence allocated to each submatrix column of the information part H i{ using 
the calculated value can be defined as 



D = {15, 15, 15, 5, 5, 5, 4, 3, 3, 3, 3, 3, 3, 3, 3} 



In the present invention, the information part H u of the parity check 
matrix H described in section A where the array code structure is described can 
be defined based on the array code structure according to the degree sequence 
defined by Equation (18). 

FIG. 7 is a diagram illustrating an example of an information part H d of 
a parity check matrix generated in the method described in section A, wherein a 
maximum variable node degree is 15. In the matrix construction illustrated in 
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FIG 7, columns and rows represent submatrix columns and submatrix rows, 
respectively. Numerals indicting the columns and rows represent cyclic shift of 
the submatrix columns and submatrix rows. Each of the numerals is comprised 
of a pxp submatrix in which p = 89. Therefore, a matrix of the information part 
5 H d generated in the forgoing method becomes a (d v p)x(d v -p) matrix. 

(II) Lifting and Hp Construction 

Generally, lifting of a matrix refers to a method of lifting a size of a 
basic matrix by submatrix replacement for a matrix having 0 and 1 in a 
10 particular position. This will be described with reference to FIG. 8. FIG. 8 is a 
diagram illustrating a 12x12 matrix obtained by matrix-lifting a basic 4x4 matrix 
in a method of replacing each element in a 4x4 matrix with a 3x3 identity matrix 
or a 3x3 0 matrix. 

1 5 As illustrated in FIG. 8, matrix lifting forms a square matrix having as 

many elements of 0 as a multiple of a target lifting size, for the elements of 0. 
Only a first row will be described. Elements of respective columns located in the 
first row becomes {0,0,1,0}, and reference numerals 801, 802, 803 and 804 are 
allocated thereto, respectively. Of all the elements, 3 elements 801, 802 and 804 

20 having a value of 0 are lifted to 3x3 matrixes 801a, 802a and 804a, respectively. 
The element 803 having a value of 1 in the first row is lifted to a 3x3 identity 
matrix 803a. The same matrix lifting is also applied to elements in other rows. 

As described above, matrix lifting generally refers to a method of 
25 extending a size of a matrix by inserting a kxk submatrix into a position of each 
element of a basic matrix composed of 0 and 1. Generally, a cyclic permutation 
matrix obtained by cyclic-shifting each column of an identity matrix is used as 
the inserted kxk matrix. 

30 A parity part H p of the parity check matrix H to be defined in the present 

invention is a generalized dual-diagonal matrix described in section B. In 
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addition, considering that H d generated in the method described in (I) of section 
C is comprised of pxp submatrixes, the H p is also constructed by lifting an rxr 
generalized dual-diagonal matrix using a pxp submatrix. 

FIG. 9 is a diagram illustrating a parity matrix H p constructed through 
matrix lifting by a pxp cyclic permutation submatrix. With reference to FIG. 9, a 
description will now be made of parity matrix lifting. A matrix lifted by the 
cyclic permutation submatrix removes 1 existing in a first row of a Jl for linear 
time encoding of a low-density parity check code and linear independency 
between respective rows. The parity part H p constructed in the manner of FIG. 9 
becomes an (rp)x(rp) matrix, and in each submatrix a h , jj denotes an offset 
value for cyclic column shift allocated to each submatrix. A description will now 
be made of a method for enabling encoding by linear calculation. 

It can be noted that when the description made in section B (Generalized 
Dual-Diagonal Matrix) is applied to the parity part H p lifted in the method of 
FIG. 9, each submatrix column of the parity part H p is selected only once by r 
calculations when 'r' is relatively prime to T . In this one selection process, it is 
possible to perform an operation of calculating a parity symbol corresponding to 
a particular column in a submatrix column. Therefore, in order to calculate parity 
symbols corresponding to all columns in a particular pxp submatrix a Jl , 
calculation for parity symbols corresponding to p columns should be performed 
only once during rp calculations. A description will now be made of a condition 
for an offset value of a submatrix for satisfying the foregoing condition. 

In the matrix construction of the parity part H p illustrated in FIG. 9, a 
relation between a row index y (,) i and a column index x (1) j in a submatrix with a 
submatrix row index = i in a diagonal line is defined as 



??=W-hi (mod/*) 



(20) 
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A relation between a row index y (2) i and a column index x (2) j in a 
submatrix with a submatrix row index = i in an offset diagonal line is defined as 



y^-^-Aw (mod/* (21) 



Based on Equation (20) and Equation (21), a value of a parity symbol 
P. can be obtained from a first row in submatrixes cr Jo and a h with a 

Jo 

submatrix row index 0 in accordance with Equation (22). 



Pk V ° (22) 



In Equation (22), v 0 denotes a partial check-sum value obtained from a 
first row of the information part H d in a parity check matrix. For a column index 
15 = j 0j a row index in a submatrix of an offset diagonal line sharing the same 
column can be expressed as 



~" A 0 1) "~/2(r- /)+! —/o /2(r-/)+l 

(23) 



20 A column index in a submatrix of a diagonal line sharing the same row 

index can be expressed as 



*(r-/) ~ fir-f) + A(r-/) " h 7»<r-/)^D + /*('-/) (24) 



25 If V and T of Equation (23) and Equation (24) are relatively prime to 

each other, all submatrixes in the parity part H p can be selected only once while 
such a process is repeated r times. Therefore, a row index in the submatrix 
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obtained again by repeating the foregoing process r times can be expressed as 

(25) 

In Equation (25), if lifting is made with all other column indexes of the 
submatrix, a variation x in column indexes updated for a submatrix a Jl existing 
in a diagonal line through r calculations is determined by 

(26) 

Theorem 2: All non-0 elements in a finite field F p = {0,1 ,2,- -,p-3 J 
defined by a particular prime number p can always be generated only once by 
performing less than 'p' additions thereon. If there is 'k' less than 'p' satisfying 
Equation (27) for a particular non-0 element 'a' in the finite field F,, in order to 
prove this, the assumption that 'p' is a prime number is disobeyed. Therefore, a 
minimum value of k satisfying Equation (27) is always 'p'. 

ka - a+ a+-- a - 0 (mod/?) 

* (27) 

Therefore, all non-0 element 'a' in the finite field F p can be generated by 
performing addition thereon less then 'p' times. From Theorem 2, if a column 
index variation of a particular submatrix o h updated by r linear calculations is 
not 0, when the j linear calculations are repeated 'p' times altogether, parity 
symbols for p column indexes can be generated. In this manner, all of rp parity 
symbols can be generated. That is, a condition of Equation (28) is satisfied. 
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(28) 

It is possible to general parity symbols corresponding to all columns in 
the submatrix by repeating r calculations of Equation (28) p times. In addition, a 
5 parity symbol corresponding to a particular column in each of (r-1) submatrixes 
existing on a diagonal line can be generated one by one during the r calculations. 
Therefore, it can be noted that parity symbols having ip different column 
indexes can be generated only once by rp linear calculations. 

10 A process of generating a matrix of a parity part H p in the method 

described above will be described with reference to FIG. 10. FIG. 10 is a 
flowchart illustrating a method for generating a matrix of a parity part according 
to a preferred embodiment of the present invention. 

15 In step 1000, a parity symbol calculation index n is set to 0 in order to 

generate a matrix of a parity part. Thereafter, in step 1002, an information 
symbol sum v 0 is calculated according to a matrix construction of an information 
part H d for a first column in a submatrix with a submatrix column index = 0, and 

a parity symbol P io for a first row in the corresponding submatrix is set to v 0 . In 

20 step 1004, a column index x (1) 0 for a parity symbol existing in a submatrix with a 
submatrix column index = 0 in a diagonal line of a parity part H p is initialized to 
Jo- 

After determining the parity symbol P io for a first row by calculating 

25 the sum v 0 of information symbols according to a matrix construction of the 
information part H d , and initializing the column index, the routine proceeds to 
step 1006. In step 1006, for a particular column index x (1) i, a row index y (2 V(r-o 
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ni a submatrix existing in an offset diagonal line sharing the same column is 
calculated by 

SiHr-f)—*i /2</4-(r-/))+l 

(29) 

Thereafter, in step 1008, a column index x (,) i+(r . 0 in a submatrix existing 
in a diagonal line having the same row index as a row index y (2) i+(r . 0 in an offset 
diagonal line is calculated by 



*"> = 4. i 

(30) 

In step 1010, a sum v (2) i+ ( r -o of information symbols existing in a row 
with a row index y (2) i+ ( r . n in a submatrix with a submatrix row index i+(r-f) is 
calculated. 



Thereafter, in step 1012, a parity symbol p 01 corresponding to a row 

•>.•.('-/» 

index x (1) i+ ( r -f) is calculated using Equation (31). 



O */ (31) 



Thereafter, in step 1014, the n is increased by 1. In step 1016, a 
submatrix row index i is updated as shown by 



r~ /+(r-2/) 

(32) 



In step 1018, it is determined whether the n is equal to rp. That is, it is 
determined in step 1018 whether encoding is completed. If it is determined in 
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step 1018 that the n is equal to rp, the routine ends, determining that encoding is 
completed. However, if it is determined in step 1018 that the n is not equal to rp, 
the routine returns to step 1006 and repeats its succeeding steps to continuously 
perform encoding. 

5 

In the foregoing process, it should be noted that all index calculations - 
submatrix row index and row index - mean a modulo-p calculation, and 
calculations on subscripts and superscripts in the index calculation process also 
mean a modulo-p calculation. Addition in the process of calculating parity 
10 symbols or information symbols is a modulo-2 calculation. An example of H p 
generated so as to satisfy such a condition is illustrated in FIG. 11. Values of 
numerals illustrated in FIG. 11 mean offset values of submatrixes. 

(IIP Overall Construction of H 

15 A parity check matrix H for a low-density parity check code to be 

defined in the present invention is defined by dividing a matrix into an 
information part H d and a parity part H p and concatenating them, for systematic 
encoding. Generation of the information part H d and the parity part H p is 
achieved in the method described in section (I) and section (II), and the 

20 generated information part H d and parity part H p become a (d v p)x(d v p) matrix 
and an (rp)x(rp) matrix, respectively. Herein, d v denotes a maximum variable 
node degree existing in an irregular low-density parity check code, 'p' denotes a 
dimension of a submatrix, and V denotes a dimension of a basic matrix given 
before lifting the parity part H p . Therefore, V should be equal to d v for 

25 concatenation. 

It can be easily understood that there is no length-4 cycle in the 
information part H d and the parity part H p of the matrix. Also, it can be easily 
understood that if a diagonal offset T of the parity part H p is set to a value 
30 approximating r/2, there is no length-4 cycle between a submatrix column of an 
information part H d where a column weight is less than r/2 and a particular 
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colunni of a parity part H p . However, in the case of columns having a maximum 
variable node degree as a weight in the information part H d , there is a possible 
case where there are columns of the parity part H p and a length-4 cycle. If a 
submatrix offset value of the parity part H p is appropriately selected, it is 
5 possible to prevent a length-4 cycle from existing in the overall parity check 
matrix H by simply removing it. 

2. Performance of LDPC Code 

Hitherto, a description has been made of a construction of a parity check 
matrix for defining a rate- 1/2 irregular low-density parity check code that can 
easily undergo systematic encoding by linear calculation. The parity check 
matrix defined in chapter 1 is divided into an information part and a parity part 
in order to generate parity symbols of a codeword according to given 
information symbols. A matrix of the information part is constructed such that a 
weight of each column in a matrix construction based on an array code structure 
has a degree approximating an optimum irregular degree distribution of a 
variable node, and a matrix of the parity part is constructed by lifting a 
generalized dual-diagonal matrix construction using a submatrix having a 
random offset. Herein, a variable node with a degree of 2 is always allocated to 
columns of the parity part. 

A description will now be made of performance of a low-density parity 
check code constructed by the parity check matrix defined in the foregoing 
description. For this, a description will first be made of a decoding algorithm and 
25 an experimental environment for evaluating performance of a low-density parity 
check code. 

A. Iterative Belief Propagation Decoding 

A low-density parity check code defined by an MxN parity check matrix 
30 can be expressed by a factor graph comprised of M check nodes and N variable 
nodes. Summarizing a message update process in a check node with a degree=d c 



10 



15 



20 
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and a variable node with a degree=d v in a log domain, message update in the 
check node can be expressed as Equation (33), message update in the variable 
node can be expressed as Equation (34), and update of a log likelihood ratio 
(LLR) can be expressed as Equation (35). 

5 

( d '~ l ^ 

l+ntanh(? ( 0) /2) 

tf-to 

l-f[tanh( 9l </) /2) 

\ <=» j 

where is a message from a check node m to a variable node nat the jth iteration and q\ n is 
a rearranged message incident to a check node m with q^ j) - q { £ at the Jth iteration. 

•••(33) 

In Equation (33), r 0 ^ is a value obtained in a j th iterative decoding 
10 process in a semi-iterative decoding process, and denotes a message delivered 
from a check node m to a variable node n. In addition, q 0) j denotes a message 
delivered from a variable node i to a check node m in a j th iterative decoding 
process. Here, T denotes a value obtained by rearranging variable nodes 
connected to a check node m from 0 to d c -l. Therefore, for i=0, it means a 
15 variable node n. 



1=1 

where q l £ is a message from variable node n to check node m at the jth iteration and q™ is 
an initial message from channel reliability. r/ ;) is a rearranged message incident to variable 
node n with r 0 (y) =r < ^ ) at the jth iteration. 

(34) 



20 In Equation (34), q° } mn is a value obtained in a j 1 iterative decoding 

process, and denotes a message delivered from a variable node n to a check node 
m. Here, r 0) j denotes a message delivered from a check node i to a variable node 
n in a j th iterative decoding process, and V denotes a value obtained by 
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rearranging check nodes connected to a variable node n from 0 to d v -l. 
Therefore, for i=0, it means a check node m. 



where 1 { „ J) is an LLR output of coded symbol n at the jth iteration and /* 0) is an initial LLR 
output from channel reliability. 



(35) 



In Equation (35), /° +1) n denotes an LLR value defined for a variable node 
n in a j" iterative decoding process. 

A description will now be made of an iterative belief propagation 
decoding process using the processes of Equation (33) to Equation (35). FIG. 12 
is a flowchart illustrating an iterative belief propagation decoding process using 
Equation (33) to Equation (35). 

15 FIG. 12 illustrates a process of decoding a received message in a receiver. 

Therefore, an initial message of a variable node n for a received message is 
defined as channel reliability of an n th symbol of a received codeword, and this 
can be expressed as 



20 Q mn = l n (36) 



n (0) 

In Equation (36), * »»» denotes an initial value of an initially-defined 

7<°> 

variable node message, and 1 « denotes an initial LLR value for an initially- 
defined variable node. In step 1200, an initial message is defined as channel 
25 reliability of an n th symbol of a received codeword as shown in Equation (36). 
Also, an iteration counter is reset. Thereafter, in step 1202, a message of a check 
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node is updated in the method given by Equation (33). In step 1204, a variable 
node and a log likelihood ratio are updated. The message of a variable node is 
updated in the method of Equation (34), and the log likelihood ratio is updated in 
the method of Equation (35). 

After the check node, variable node, and log likelihood ratio are all 
updated in the foregoing process, a value of the updated log likelihood ratio 
undergoes hard decision in step 1206. Thereafter, in step 1208, parity check is 
performed on a received message based on the hard-decided value. If the parity 
check result has a value of 0 meaning that decoding has been successfully 
performed, the decoding ends in step 1216. However, if the parity check result 
does not have a value of 0, it is determined in step 1210 whether iteration has 
reached a predetermined iteration number. If it is determined in step 1210 that 
iteration has reached the predetermined iteration number, a decoding failure 
process is performed in step 1214, determining that although decoding is 
performed additionally, there is low probability that the decoding will be 
achieved successfully. However, if it is determined in step 1210 that iteration has 
not reached the predetermined iteration number, the iteration counter is increased 
by 1 in step 1212, and then, the routine returns to step 1202. 

B. Simulation Environment 

An experimental environment for evaluating performance of a low- 
density parity check code constructed by a parity check matrix defined in the 
present invention is illustrated in Table 1 . 

Table 1 

- Code rate = 1/2, maximum variable node degree =15 

- Offset (f) in Hp = 7 

- Frame size = 435(p=29), 795 (p=53), 1545(/?=103), 3855(p=257) 

- Binary antipodal signaling over AWGN channel 

- Iterative belief propagation decoding 
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- Floating point simulation 

- Maximum number of iterations = 160 

- Stop by parity check at each iteration 

- Frame error rate (FER) & information bit error rate (BER) evaluation 



For convenience, the present invention limits a code rate of a low- 
density parity check code to 1/2. However, parity check matrix design of a low- 
density parity check code for supporting various code rates can be discussed as a 
method for lifting a scheme proposed in the present invention according to a 
code rate. A maximum variable node degree is set to 15 in a factor graph of a 
low-density parity check code by a parity check matrix defined in the present 
invention, in order to acquire excellent performance and prevent a hardware size 
from increasing in implementation. An offset f in a matrix of a parity part H p is 
set to 7 in order to simply remove a length-4 cycle from a parity check matrix, in 
the maimer described in (III) of section C. A frame size of a low-density parity 
check code simulated in the present invention is set to a value similar to an 
encoder packet (EP) size specified in the cdma2000 lxEV-DV standard so that 
performance of the low-density parity check code is compared with performance 
of a turbo decoder currently used in the cdma2000 lxEV-DV standard. In a 
decoding process of a turbo code for performance comparison, a log-MAP 
algorithm is used, and the maximum number of iterations is limited to 8. 

The maximum number of iterations of an iterative belief propagation 
decoder for a low-density parity check code is set to 160, taking into 
consideration the fact that a convergence speed of an irregular low-density parity 
check code is relatively low in the general density evolution technique. For 
experiments on the decoder, floating point simulation expressing extrinsic & 
LLR information in real values was performed. In a check node update process, 
in order to prevent overflow of a function (p(x)=-log(tanh(|x|/2), a value of <p(x) 
is limited to 20 for x in |x| <10" 8 . A decoder for a low-density parity check code 
performs parity check at each iteration and adopts the parity check result as a 
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stopping criterion, so that no error is detected in a parity check process. However, 
when an error occurs actually, the decoder classifies the defective frame as a 
frame having an undetected error. Finally, a frame error rate and a bit error rate 
are used as measures for evaluating performance of the low-density parity check 
5 code designed in the present invention. The frame/bit error is considered only 
when an error occurs in information symbols. A description will now be made of 
performance of the parity check code simulated in the foregoing simulation 
environment. 

10 C. Simulation Results 

FIG 13A is a diagram illustrating an example of a matrix of an 
information part H d for n=870 and p=29, FIG 13B is a diagram illustrating an 
example of a matrix of a parity part H p for n=870 and p=29, and FIG 13C is a 
simulation result graph illustrating a comparison between a low-density parity 

1 5 check code and a turbo code. 

The matrix of an information part and the matrix of a parity part 
illustrated in FIGs. 13 A and 13B are constructed in the foregoing method 
according to the present invention, and actually, the two matrixes are 

20 concatenated before being encoded. A description will now be made of the 
conditions for a case where the performances illustrated in FIG 13C are 
available. For performance comparison between two coding methods, 
performance for a case where an encoder packet (EP) size of the cdma2000 
lxEV-DV standard is 408 and a code rate is 1/2 is shown together. As illustrated 

25 in FIG 13C, it can be understood that for a small block size, performance of a 
turbo code used in the cdma2000 lxEV-DV standard is slightly superior to 
performance of the low-density parity check code proposed in the present 
invention in terms of both a frame error rate (FER) and a bit error rate (BER). 
This shows that a parity check matrix construction defining the low-density 

30 parity check code defined in the present invention is not an optimized 
construction in terms of decoding performance of a codeword with a small size. 
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That is, it can be noted that when a coded-block size of a low-density parity 
check code is small, in order for the low-density parity check code to show 
better performance as compared with the existing known turbo code, an 
optimization operation on the parity check matrix should be performed further. 
5 However, it is shown that the low-density parity check code is superior to the 
turbo code in terms of decoding performance of a codeword with a large block 
size. 

Table 2 illustrates the average number of iterations for each Eb/No for 
10 n=870. According to Table 2, although a maximum of 160 iterations are 
performed on a given low-density parity check code, an actual decoding result 
shows that it is possible to obtain sufficient performance with only less than 10 
iterations considering a high SNR within 20 iterations. 

15 Table 2 



Eb/No (dB) 


Average number of iterations 


1.6 


15.085 


1.8 


11.353 


2.0 


9.317 


2.2 


8.01 


2.4 


7.087 



FIG. 14A is a diagram illustrating an example of a matrix of an 
information part H d for n=1590 and p=53, FIG. 14B is a diagram illustrating an 
example of a matrix of a parity part H p for n=1590 and p=53, and FIG 14C is a 
20 simulation result graph illustrating a comparison between a low-density parity 
check code and a turbo code. 

The matrix of an information part and the matrix of a parity part 
illustrated in FIGs. 14A and 14B are constructed in the foregoing method 
25 according to the present invention, and actually, the two matrixes are 
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concatenated before being encoded. 

A description will now be made of the conditions for a case where the 
perforaiances illustrated in FIG. 14C are available. For performance comparison 

5 between two coding methods, performance for a case where an encoder packet 
size of the cdma2000 lxEV-DV standard is 792 and a code rate is 1/2 is shown 
together. It is noted from FIG. 14C that as a block size of a low-density parity 
check code increases, the low-density parity check code is almost equal to the 
existing turbo code in frame error rate. However, the low-density parity check 

10 code is still inferior to the existing turbo code in bit error rate. This is because 
the turbo code is lower than the low-density parity check code in terms of the 
ratio of bit errors existing in a defective frame of the turbo code. That is, the 
turbo code is lower than the low-density parity check code in BER/FER for the 
same Eb/No. In this case, it means that the turbo code is lower than the low- 

15 density parity check code in average number of bit errors existing in a defective 
frame. However, in an actual wireless communication system employing Hybrid 
Automatic Retransmission Request (ARQ), because it is more preferable that 
FER performance is superior to BER performance, it can be considered from 
FIG. 14C that for n=l590, the low-density parity check code according to the 

20 present invention is equal in performance to the turbo code used in the 
cdma2000 lxEV-DV standard. 

Table 3 illustrates the average number of iterations for each Eb/No for 
n=1590. As a block size increases, the average number of iterations is slightly 
25 higher than that in Table 2. However, it can be understood that at a sufficiently 
high SNR, sufficiently high performance can be obtained with only less than 20 
iterations. 



Table 3 



Eb/No (dB) 


Average number of iterations 


1.1 


40.308 
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1.3 


21.504 


1.5 


14.858 


1.7 


12.018 


1.8 


11.020 



FIG. 15A is a diagram illustrating an example of a matrix of an 
information part H u for n=3090 and p=103, FIG. 15B is a diagram illustrating an 
example of a matrix of a parity part H p for n=3090 and p=103, and FIG. I5C is a 
simulation result graph illustrating a comparison between a low-density parity 
check code and a turbo code. 

The matrix of an information part and the matrix of a parity part 
illustrated in FIGs. 15A and 15B are constructed in the foregoing method 
according to the present invention, and actually, the two matrixes are 
concatenated before being encoded. 



A description will now be made of the conditions for a case where the 
performances illustrated in FIG. 15C are available. For performance comparison 
between two coding methods, performance for a case where an encode, packet 
size of the cdma2000 lxEV-DV standard is 1560 and a code rate is 1/2 is shown 
together. It is noted from FIG. 15C that as a block size of a low-density parity 
check code increases, the low-density parity check code is superior to the 
existing turbo code with a similar size in frame error rate. Although the low- 
density parity check code is still inferior to the existing turbo code in bit error 
rate as illustrated in FIG. 14C, a difference between them is greatly reduced, 
closely approximating performance of the turbo code. 

Table 4 illustrates the average number of iterations for each Eb/No for 
n=3090. As a block size increases, the average number of iterations is slightly 
higher than that in Table 2. However, it can be understood that at a sufficiently 
high SNR, sufficiently high performance can be obtained with only less than 20 
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iterations. 



Table 4 



Eb/No (dB) 


Average number of iterations 


1.1 


38.495 


1.2 


22.051 


1.4 


16.190 


1.5 


14.623 



FIG. 16 A is a diagram illustrating an example of a matrix of an 
information part H d for n=7710 and p=257, FIG. 16B is a diagram illustrating an 
example of a matrix of a parity part H p for n=7710 and p=257, FIG. 16C is a 
simulation result graph illustrating a comparison between a low-density parity 
check code and a turbo code, and FIG 16D is a simulation result graph for a 
low-density parity check code according to a variation in the number of 
iterations. 

The matrix of an information part and the matrix of a parity part 
illustrated in FIGs. 16A and 16B are constructed in the foregoing method 
according to the present invention, and actually, the two matrixes are 
concatenated before being encoded. 

A description will noAv be made of the conditions for a case where the 
performances illustrated in FIG. 16C are available. For performance comparison, 
performance for a case where an encoder packet size of the cdma2000 lxEV-DV 
standard is 3864 and a code rate is 1/2 is shown together. It is noted from FIG 
16C that when a block size of a low-density parity check code is very large, the 
low-density parity check code is very superior to the turbo code with a similar 
size in performance. In particular, it can be understood that in the case of a turbo 
code with a frame length of 3864 used in the cdma2000 lxEV-DV standard, both 
a frame error rate and a bit error rate have an error floor at a high SNR. However, 
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it is observed that in the case of the low-density parity check code defined in the 
present invention, such an error floor does not occur. Such performance is very 
advantageous to a communication system employing H-ARQ. Accordingly, it 
can be noted that the low-density parity check code defined in the present 
invention is very superior in performance to the turbo code used in the existing 
standard. Table 5 illustrates the average number of iterations for each Eb/No for 
n=7710. 



Table 5 



Eb/No (dB) 


Average number of iterations 


0.9 


41.036 


1.0 


31.034 


1.1 


25.507 


1.2 


22.170 



In Table 5, as a block size of the low-density parity check code increases 
the average number of iterations increases at high rate as compared with that of 
Table 4. That is, the low-density parity check code supporting a maximum frame 
size similar to the size defined in the cdma2000 lxEV-DV standard requires 20 
or more iterative decodings even at a high SNR. Therefore, in order to show 
optimal performance with minimum delay, it is necessary to observe 
performance of an LDPC decoder according to the number of iterations. 

FIG. 16C is a diagram illustrating available FER/BER performance of 
the low-density parity check code having the constructions shown in FIGs. 16A 
and 16B in which the maximum number of iterations is limited to 40, 80, 120 
and 160. Unlike this, if the maximum number of iterations is limited to 40, 80 
and 120, slight performance deterioration occurs. However, except the case 
where the maximum number of iterations is limited to 40, the performance 
deterioration is not considerable. Therefore, although the maximum number of 
iterations is limited to about 80, the performance deterioration problem is not 
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considerable. In particular, when a block size is smaller than this, the 
performance deterioration caused by limiting the maximum number of iterations 
to 80 is further reduced, so that it becomes ignorable. Therefore, the low-density 
parity check code proposed in the present invention can obtain sufficient 
5 performance during decoding, if the maximum number of iterations is set to 
about 80, on the assumption that the maximum block size is 7710. 

3. Conclusion 

The present invention defines a low-density parity check code capable of 
10 enabling efficient encoding and obtaining excellent decoding performance only 
for a code rate = 1/2. To this end, the present invention partitions a parity check 
matrix defining a low-density parity check code into two parts, and defines a 
part corresponding to information on a codeword in the parity check matrix 
using an array code structure whose column weight is higher than 2. In addition, 
1 5 the present invention defines a part corresponding to a parity of a codeword in 
the parity check matrix as a generalized dual-diagonal matrix whose column 
weights are all 2. The parity check matrix defined in this manner generates an 
irregular low-density parity check code whose maximum variable node degree 
d v =15, in which a threshold for error free is 0.9352 in an AWGN channel by the 
20 density evolution technique based on Gaussian approximation. This shows 
performance approximating Shannon channel capacity by 0.5819 dB for a rate- 
1/2 code. 

Based on the fact that linear encoding is available, once a parity check 
25 matrix is designed in a random scheme without considering possibility of actual 
linear encoding, an excellent low-density parity check code approximating the 
Shannon threshold by 0.5819 dB or closer can be generated. In this case, the 
encoding process becomes complicated, making it difficult to enable actual 
implementation. Because enabling linear encoding serves as one restrictive 
30 condition for a parity check matrix defining a low-density parity check code, the 
parity check matrix should have a random structure. Therefore, a possibility for 
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showing excellent performance will be lower than that of a parity check matrix 
that is randomly defined without considering linear encoding. Thus, application 
of the present invention shows relatively excellent performance considering the 
linear encoding. That is, the present invention shows relatively excellent 
5 performance despite a restriction on the parity check matrix due to the linear 
encoding. 

In order to show excellent decoding performance, a size of a codeword 
is increased and a maximum variable node degree is increased accordingly in 

10 defining a degree distribution of a parity check matrix, making it possible to 
design a low-density parity check code having performance approximating the 
Shannon capacity limit. However, in this case, because an increase in the 
maximum variable node degree increases complexity of the decoder during 
decoding, it cannot contribute to actual implementation. Therefore, that decoder 

15 complexity is allowable means that a reasonable maximum variable node degree 
is set in implementation. 

The low-density parity check code by a parity check matrix defined in 
the present invention has a generalized dual-diagonal matrix in which a parity 

20 part of a codeword in the parity check matrix is lifted by cyclic permutation of a 
pxp identity matrix. Therefore, it can be simply encoded by simple linear 
calculations. In addition, when decoding performance by iterative belief 
propagation was compared with decoding performance of the turbo code used in 
the cdma2000 lxEV-DV standard, it was proven through simulations that a 

25 lower frame error rate can be obtained for the similar frame size. In particular, 
compared with the turbo code, the low-density parity check code according to 
the present invention has a lower frame error rate, and when it is used together 
with the Hybrid -ARQ technique, it will serve as a point of excellence. 
Therefore, the low-density parity check code defined in the present invention can 

30 be simply encoded by simple linear calculations. Also, in implementation of the 
decoder, check node processors and variable node processors are implemented in 
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parallel in rows and columns of respective submatrixes, thereby enabling fast 
decoding. 

In the low-density parity check code defined in the present invention, it 
5 can be simply checked that there is no length-4 cycle in a factor graph defining 
the low-density parity check code according to a characteristic of arrangement of 
the submatrixes constituting a parity check matrix. 

As described above, application of the present invention can provide 
10 performance similar to or superior to that of the turbo decoder, and in particular, 
can generate a low-density parity check code capable of reducing a frame error 
rate. In addition, the low-density parity check code according to the present 
invention can reduce decoding complexity. 



WO 2005/036758 



PCT/KR2004/00263H 



-41 - 



WHAT IS CLAIMED IS: 

1. A method for generating a low-density parity check code 
consisting of an information-part matrix and a parity-part matrix, comprising the 
steps of: 

changing the information-part matrix to an array code structure, and 
allocating a degree sequence to each of submatrix columns; 

extending the parity-part matrix such that an offset value between 
diagonal lines has a predetermined value in a generalized dual-diagonal matrix 
which is the parity-part matrix; 

lifting the generalized dual-diagonal matrix; 

determining an offset value for cyclic column shift for each submatrix of 
the lifted generalized dual-diagonal matrix; and 

performing an encoding process for determining a parity symbol 
corresponding to a column of the parity-part matrix. 

2. The method of claim 1, wherein the degree sequence is formed 
in accordance with Equation (37). 

D = {15, 15, 15, 5, 5, 5, 4, 3, 3, 3, 3, 3, 3, 3, 3} 

(37) 

3. The method of claim 1, wherein the offset value between 
diagonal lines is relatively prime to the number of columns in the generalized 
dual-diagonal matrix. 

4. The method of claim 1, wherein the number of rows in the 
submatrix is a prime number. 

5. The method of claim 1, wherein a difference between a sum of 
offset values for cyclic row shift of a submatrix on a diagonal line in a 
generalized dual-diagonal matrix which is the parity-part matrix and a sum of 
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offset values for cyclic row shift of a submatrix on an offset diagonal line is not 
0. 

6. The method of claim 1, wherein the encoding process comprises 
5 the process of: 

(a) determining a parity symbol of a first row in a submatrix with a 
submatrix column index 0 on a diagonal line of the parity-part matrix; 

(b) setting a row index in a submatrix of a parity symbol being identical 
to the determined parity symbol in column index in a submatrix in a submatrix 

10 on an offset diagonal line having the same submatrix column index as the 
submatrix column index of the set parity symbol; 

(c) determining a parity symbol having the same row index in the set 
submatrix in a submatrix on a diagonal line having the same submatrix row 
index as the submatrix row index of the submatrix on the offset diagonal line; 

15 and 

(d) repeatedly performing the steps (b) and (c) until generation of the 
parity matrix is completed. 

7. The method of claim 6, wherein in step (a), the parity symbol is 
20 determined by a sum of information symbols of the information-part matrix 

existing in the same row as a row index in the submatrix whose parity symbols 
are determined. 

8. The method of claim 6, wherein in step (b), the row index in the 
25 submatrix is set in accordance with Equation (38). 

(38) 

where v (2) i+(r . {) denotes a row index in a submatrix with a submatrix column 
index i on an offset diagonal line, x (1) j denotes a column index in a submatrix 
with a column index i existing in a diagonal line, and j2(i+(r-f))+i denotes an offset 
30 value for cyclic column shift of a submatrix with a submatrix column index i on 
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the offset diagonal line. 

9. The method of claim 6, wherein in step (c), the parity symbol is 
determined in accordance with Equation (39). 

•\-w) *i Yh*-i\ (39) 

where PjLXX denotes a parity symbol corresponding to x (1) i+(r . n? p tl) denotes a 
parity symbol corresponding to a column index x (1) j, and v (2) denotes a sum 

of information symbols existing in a row with a row index y (2) i+( r . 0 in a submatrix 
with a submatrix column index i+(r-f). 
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Example of Hj with irregular distribution of d v = 15 (p=89) 
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Example of Hj with irregular distribution of d v = 15 (p=89) 
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